Dynamic vortex mass in clean Fermi superfluids and superconductors 
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We calculate the dynamic vortex mass for clean Fermi superfluids including both s- and d-wave 
superconductors as a response to a vortex acceleration. Assuming a finite quasiparticle mean free 
time, the vortex mass appears to be a tensor. The diagonal component dominates in the limit of 
long mean free time while the off-diagonal mass takes over in the moderately clean regime. 



The vortex mass in superfluids and superconductors is 
a long standing problem in vortex physics and remains 
to be an issue of controversies. There are different ap- 
proaches to its definition. One approach consists in cal- 
culating the vortex free energy increase due to the vortex 
velocity 0. First used by Suhl §] this approach yields 
the mass of the order of one quasiparticle mass (elec- 
tron, in case of superconductor) per atomic layer (see 
Rcf. for an extensive review). Another approach is 
based on finding the force necessary to support an un- 
steady vortex motion. Identifying then the contribution 
to the force proportional to the vortex acceleration, one 
defines the vortex mass as a coefficient of proportional- 
ity. This method was first applied for vortices in super- 
clean superconductors in Ref. B and since then was used 
widely (see for example, |||-|7|]). The resulting mass is of 
the order of the total mass of all particles within the area 
occupied by the vortex core. We will refer to this mass 
as to the dynamic mass. The dynamic mass originates 
from the inertia of excitations localized in the vortex core 
and can also be calculated as the momentum carried by 
localized excitations ||. It is much larger and thus it is 
much more important than the mass obtained from the 
energy considerations. 

In the present Letter we develop a regular microscopic 
approach for calculating the dynamic vortex mass in a 
general case of a finite relaxation time r of nonequilib- 
rium excitations produced by the moving vortex. Using 
the Boltzmann kinetic equation for quasiparticles local- 
ized in the vortex core, we derive the equation for the 
vortex dynamics which contains the inertia term together 
with all the forces acting on a moving vortex. This ap- 
proach is applied to both s- and d-wave superconductors. 
We find that dynamic mass displays a novel feature: it 
is a tensor whose components depend on the quasipar- 
ticle relaxation time. In s-wave superconductors, this 
tensor is diagonal in the superclean limit. The diagonal 
mass decreases rapidly as a function of the mean free 
time, and the off-diagonal components dominate in the 
moderately clean regime. In d-wave superconductors, the 
transition from a diagonal mass (in the extremely clean, 
high-r limit) to an off-diagonal mass tensor (in the mod- 
erately clean, low-r regime) occurs via the intermediate 



universal regime of the flux flow predicted in || (see also 
[|l0|,[ri]] ) , where both components are of the same order 
of magnitude. Our results agree with the previous work 
in the limit r -> oo. 
Forces on a moving vortex. The Boltzmann kinetic 
equation for quasiparticles localized in the vortex core is 
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Here e n (a,/j,) is the quasiclassical spectrum of a quasi- 
particle, which depends on canonically conjugated vari- 
ables a and [i where a is the angle of the particle veloc- 
ity vj_ = ±p_i_ in the plane perpendicular to the vortex 
axis, and \x = z fp±b is the "angular momentum" defined 
through the impact parameter 6 of a particle. The impact 
parameter is a good quantum number in the quasiclassi- 
cal limit and so is fi, even if there is no axial symmetry. 
The =F signs refer to particles and holes in the sense of 
the normal-state Fermi surface with positive and nega- 
tive effective masses (in the plane perpendicular to the 
vortex axis), respectively. We take the coordinate z axis 
along the vortex circulation z = sgn(e)B/£? where B is 
the magnetic induction. It can be shown |12] that, for the 
vortex dynamics, the Boltzmann-equation approach |13f| 
is equivalent to the full-scale Green function formalism 
[j]~^ , P~5|] . Eq. (|l]) can be derived for slow in time variations 
of the distribution function such that the characteristic 
w is small as compared to T or to the bulk gap Ao. 

For simplicity, we restrict ourselves to the case of well 
separated vortices, H <C H C 2, and consider the limit 
lo c t <C 1 where lo c = \e\H/m*c is the cyclotron frequency. 
The limit lo c t <C 1 is the most realistic for superconduc- 
tors; it also corresponds to uncharged superfluids. In this 
limit, the delocalized quasiparticles above the gap are in 
equilibrium with the heat bath and do not participate in 
the vortex motion . 

We seek the distribution function in the form / = 
/ (0) (e«) + fi, where /(°)(e) = 1 - 2n(e) and n(e) is 
the Fermi function while f\ describes the deviation from 
equilibrium. For a vortex moving with the velocity 
u the angular momentum is time dependent: fi(t) = 
[(r — ut) x p] • z. The driving term in Eq. (Q) is 
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The spectrum of bound states in the vortex core has an 
anomalous chiral branch [|l6| eo(/Lt) whose energy varies 
from — Ao to Ao crossing zero as the impact parame- 
ter of a particle in the core varies from — oo to +oo. 
The anomalous branch corresponds to the radial quan- 
tum number n = 0. The energies of states with n ^ 
are separated from e — by a gap of the order of Ao. 
Defining uj n = =F<9e„/9/i, we ensure u n to be positive on 
the branch n — for both particles with vj_ = p^ and 
holes with vj_ = — pj_. Note that it is the Fermi veloc- 
ity v_l rather than the momentum pj_ that enters the 
quasiclassical equations which determine the excitation 
spectrum. Consequently, it is de n /db that has the same 
sign for particles and holes. 

Multiplying Eq. (Q) by p±/2 and summming up over 
all the quantum numbers, we obtain 
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where the l.h.s. of Eq. (|3|) is the force from the envi- 
ronment on a moving vortex derived in [fj"5[ . It is the 
momentum transferred from excitations to the vortex: 
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The first term in the r.h.s. of Eq. (g) is the force 
exerted on the vortex by the heat bath via excitations 
localized in the vortex core: 
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Its component parallel to u gives the friction force. The 
transverse component is nothing but the spectral flow 
force due to the localized excitations jL7|. The second 
term in the r.h.s. of (||) is the change in the vortex mo- 
mentum defined as 
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Finally, the third term in the r.h.s. of Eq. (||) is ob- 
tained from Eq. (j^) (we assume here an s-wave pairing 
for simplicity). The sum over n contains only the spec- 
tral branch with n = because de n /d/i is odd in /x for 
n =/= 0. It can be written as 
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Here we put N = N e -N h and N s = N — N n where N n is 
the normal density. The term containing N n in Eq. (0) 



is the Iordanskii force JL8| while the third one is a part 
of the spectral flow force due to the excitations above 
the gap |l^]. The last two terms in Eq. (0) together 
can be written as the transverse force produced by the 
delocalized excitations scattered by the vortex potential 
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with the transverse cross section obtained in Ref. pi 
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Note that the similar longitudinal force can be neglected 
since it is smaller than the transverse one by a factor 
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The equation of the vortex dynamics is obtained by the 
variation of the superfluid free energy plus the energy in 
the external field with respect to the vortex displacement. 
The variation of the superfluid free energy gives the force 
from the environment, F onv , while the variation of the en- 
ergy in the external field produces the external force. For 
charged superfluids the external force is just the Lorentz 
force @F L = ($o/c)[j s x z]sgn(e) where $ = 7rc/2|e| 
is the magnetic flux quantum. This formula holds also 
for uncharged systems since the charge drops out giving 
the external force in the form Fl = "kN s \v s x z]. In 
the absence of pinning, Fl + F cnv = since the total 
energy is translationally invariant. Using Eq. (JlTJ) the 
force balance takes the form 
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In a Galilean invariant system, where N is the to- 
tal number of carriers, the 7V s -term can be combined 
with the Lorentz force to give the Magnus force Fm = 
7riV s [(v s — u) x z]. Eq. ( |Tl| ) becomes 
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The physical meaning of the Eq. ([lT]) [or (fL2|) l is sim- 
ple. The equation accounts for all the forces acting on a 
moving straight vortex line: the Magnus force from the 
superfluid component, the force due to the scattering of 
normal excitations F sc , and the force from the heat bath 
transferred to the vortex by localized excitations, Fi oc . 
The r.h.s. of Eqs. ( pT[ |l2] ) comes from the inertia of 
excitations localized in the vortex core and is identified 
as the change in the vortex momentum. The definition 
of Eq. ^ is similar to that used in Refs. HQ. Note 
that the delocalized quasiparticles do not contribute to 
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the vortex momentum in the limit lu c t <C 1 because their 
distribution is equilibrium ]l5[ | with fx = 0. 

Vortex mass. Having defined the vortex momentum, 
we calculate the vortex mass. The distribution function 
of localized particles is M| 



h = - 
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where 7h and 70 are to be found from Eq. (Q). The 
vortex momentum becomes Pi = MikUk] it has both lon- 
gitudinal and transverse components with respect to the 
vortex velocity. 

For a vortex with the symmetry not less than the four- 
fold, the effective mass tensor per unit length is = 
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and the same expression for Mj_ e ^h with 7h replaced with 
70 • The e, h subscripts indicate the integrations over the 
electron and hole parts of the Fermi surface, respectively. 
Only the branch with n = gives the contribution to the 
transverse mass because 70 is odd in /1 for n^O. 

Eq. (^|) does not contain any explicit contribution from 
the so called "backflow mass" . The backflow vortex mass 
was introduced pcfl to account for an apparent similar- 
ity between a vortex core and an impenetrable cylinder 
immersed into a fluid. A hydrodynamic dipole counter- 
flow would arise around the cylinder moving with respect 
to the fluid giving rise to the associated backflow mass 
EjJ . We see, however, that the total inertial force comes 
from the change in the momentum of localized excita- 
tions only. In a sense, Eq. (||) takes care of all sources of 
the vortex mass including the possible backflow effect. 

We also comment on arguments of Ref. Q that the 
Coulomb interactions from charge variations produced 
by the moving vortex would screen the dynamic mass in 
charged superfluids. It is well known that moving vor- 
tex induces a non-equilibrium chemical potential which 
compensates the changes in the charge density (see, for 
example, ^,^3|) . Therefore, the vortex motion does 
not violate the charge neutrality and the abovementioncd 
screening effects are absent. (Note also that such a 
screening would first of all have screened out the whole 
vortex response to an electric current.) 

s-wave superconductors. For s wave superconductors, 
uj n is independent of a. If the collision integral is taken in 
the T-approximation (<9//9i) coll — — /i/r n , the solution 
to Eq. (|l|) for a steady vortex motion is [^5) 
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where uj n ~ T 2 /Ep. At low temperatures one has 
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where (...) = (ir/Vp) J p\ dp z (. . .) is the averaging over 
the Fermi surface with the volume Vp. The mass is of 
the order of the mass of electrons in the area occupied 
by the vortex core M ~ Tr^mN. In the superclean limit 
T 2 r/Ep 3> 1 where uj n T ^> 1 the mass tensor is diagonal 
Mik — Mii^jfe, and we obtain the known result @,|||§: 
M|| e , h = TrN e . h (u- 1 ). 

The mass decreases with r. In the moderately clean 
regime T 2 r/Ep 1 where oj n T <§; 1, the diagonal com- 
ponent vanishes as r 2 , and the mass tensor is dominated 
by the off-diagonal part. If tq is independent of energy 
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The vortex mass is not a constant quantity for a given 
system: it may depend on the frequency u> of the ex- 
ternal drive. Indeed, for lu comparable with uj q and r — 1 , 
both the forces and the effective mass acquire a frequency 
dispersion. From Eq. (fil) one finds 
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w « T n + (1 - «WT„) 2 ' 70 W 2 T 2 + (1 - iuJT n ) 2 ' 

d-wave superconductors. We take the d-wave order 
parameter in the form A = A (/9, 4>) sin (2a) e^. The 
interlevel spacing w„(a) for the quasiclassical spectrum 
e n (a,/i) of particles in the vortex core depends now on 
the momentum direction a. 

One distinguishes three different regimes for systems 
with different purity. In the order of increasing r we 
have (i) the moderately clean regime, (ii) the universal 
regime, and (hi) the extremely clean limit. The moder- 
ately clean limit corresponds to T 2 r/Ep <C 1. In this 
case, the factors 7 in Eq. ( |l3| ) are 7h = oj 2 (a)r 2 and 
70 = uj n (a)T n . The mass tensor has qualitatively the 
same behavior as for an s-wave superconductor. 

In the superclean limit T 2 r/Ep 3> 1, the specifics of a 
d-wave superconductor is more pronounced PJTl|. Con- 
sider temperatures T <C T C ^J H/ H C 2- At such tempera- 
tures, only localized states on the anomalous branch are 
excited. We need the quasiclassical spectrum for momen- 
tum directions near the gap nodes. The distance between 
the quasiclassical levels on the anomalous branch near the 
gap nodes is [O 
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where Q.q ~ T 2 /Ep, and oj c = \e\Hv±/cp± is the cy- 
clotron frequency. For comparison, lu c /CIq ~ H/H C 2- 
The factors 7h and 70 in Eq. (|l3j ) satisfy the set of 
equations which follow from Eq. (|l[) [9|]Tl|] 
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Here U(a) = [wo(a)ro] . The masses are 
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where F (a) = J " U(a') da' and r = ro(a = 0). The in- 
tegral over da is determined by the angles a ~ H / H C 2- 
For this range, the terms with U(a) in Eq. Jl7| ) dominate 
and the factors 7h, o are 
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Here 
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and Eq = y/flouj c is the distance between the exact quan- 
tum levels in the vortex core. They can be found by 
semi-classical quantization of the azimuthal motion of 
the particle having the spectrum eo(a, fx) |0]. Finally, 
the masses become 
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The universal regime |9j corresponds to A > 1, i.e., 
Eqt <C 1. In this limit, the spacing between the ex- 
act quantum levels is much smaller than the relaxation 
rate, while the average distance between the quasiclas- 
sical levels is larger than the relaxation rate, l^r 3> 1. 
Here, the both masses are of the same order of magni- 
tude Mil ~ M± ~ Nr. In the extremely clean limit when 



II 

A <C 1, i.e., E t 3> 1, one obtains Mm e j } 
which agrees with the result of Ref . ||pT| . 



= TrN e<h (E^ 1 ) 
The transverse 
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In conclusion, we have calculated the dynamic vortex- 
mass tensor as a response to a slow vortex acceleration. 
The diagonal component dominates in the limit of a very 
long mean free time. In the opposite limit of fast re- 
laxation, the off-diagonal mass is more important. The 
order of magnitude of the dynamic mass is the mass of 
all the quasiparticles in the area occupied by the vortex 
core. This mass is much larger than both the electro- 
dynamic and the Suhl masses and is associated with the 
inertia of excitations localized in the vortex core. 
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